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ABSTRACT 

SONYC - Suhstellar Objects in Nearby Young Clusters ~ is a survey program to investigate the 
frequency and properties of suhstellar ohjects with masses down to a few times that of Jupiter in 
nearby star- forming regions. For the ~lMyr old p Ophiuchi cluster, in our earlier paper we reported 
deep, wide-field optical and near-infrared imaging using Subaru, combined with 2MASS and Spitzer 
photometry, as well as follow-up spectroscopy confirming three likely cluster members, including a new 
brown dwarf with a mass close to the deuterium-burning limit. Here we present the results of extensive 
new spectroscopy targeting a total of '-^lOO candidates in p Oph, with FMOS at the Subaru Telescope 
and SINFONI at the ESO's Very Large Telescope. We identify 19 objects with effective temperatures 
at or below 3200 K, 8 of which are newly identified very-low-mass probable members of p Oph. Among 
these eight, six objects have Tog < 3000K, confirming their likely suhstellar nature. These six new 
brown dwarfs comprise one fifth of the known suhstellar population in p Oph. We estimate that the 
number of missing suhstellar objects in our survey area is ^ 15, down to 0.003 — 0.03711© and for 
Ay = 0—15. The upper limit on the low-mass star to brown dwarf ratio in p Oph is 5.1 ± 1.4, while 
the disk fractions are ~ 40% and ~ 60% for stars and BDs, respectively. Both results are in line with 
those for other nearby star forming regions. 

Subject headings: stars: circumstellar matter, formation, low-mass, brown dwarfs - planetary systems 



1. INTRODUCTION 

Understanding the origin of the stellar initial mass 
function (IMF) is one of the major topics in astrophysics. 
The low-mass end of the IMF, in particular, has been 
subject of numerous observ ational and theoreti cal stud- 
ies over the past decade (see Bonnell et al.|2067 1 . Brown 
dwarfs (BD; M < 0.08 M©) are an important class of ob- 
jects to test the mass dependence in the formation and 
early evolution of stars and planets. Establishing a cen- 
sus of young suhstellar objects is a fundamental prereq- 
uisite to test the competing theoretical scenarios trying 
to explain their origin. These scenari os iii clude turbulent 
fragmentation (Pad oan & Nordlund 20Q4|), ejections from 
multiple systems ( Bate 2009) or proioplanetary disks 
(Stamatellos & Whitworth 2009). Further issues that 
can be addressed by investigating young brown dwarfs 
are, among others, the evolution of accretion disks (e.g., 
the impact o f photoevap oration) , the formation of multi- 
ple systems ( Bate||2009[), the mass dependence of planet 
formation (Payne & Lodato 2007), and the early reg- 
ulation of angular momentum ( [Scholz fc Eisl ofFcf 2004, 
2005 ) . Young BDs also provide exc ellent test- beds for 
low-g ravity atmosphere models (e.g. Kirkpatrick et al. 

2om. 

Surveys in star-forming regions have revealed the ex- 
istence of free-floating objects with masses below the 
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Deuterium-burning limit (^ . 015 M©; Zapatero Osorio 
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2000[ [Lucas et aL |20051 [Jayawardhana fe Ivanov 



Luhman et al. 2008a), Te^ comparable to tlie 



masses ot giant planets! Recently, Sumi et al. (2011) 
presented evidence for the existence of a population of 
unbound or distant Jupiter-mass objects in the field, 
that may have formed in proto-planetary disks and sub- 
sequently scattered. Nevertheless, current surveys in star 
forming regions are incomplete in the suhstellar regime, 
which is especially pronounced in the planetary-mass do- 
main. Some of the deepest surveys are biased by design, 
either because they are based on mid-infrared data and 
thus will only find objects surrounded by dusty disks 



(e.g., Luhman et al. 2008a) or because they are carried 



Haisch et al. 



out using narrow- band methane imaging, which probes 
only a specific range of spectral types (e.g., " 
20T0| ). 

^DNYC - Suhstellar Objects in Nearby Young Clus- 
ters - is an ongoing project to provide a complete census 
of the brown dwarf and planetary-mass object (PMO) 
population in nearby young clusters and to establish the 
frequency of suhstellar mass objects as a function of clus- 
ter environment. The survey is based on our own ex- 
tremely deep optical- and near-infrared imaging, com- 
bined with the Two Micron All Sky Survey (2MASS) and 
Spitzer photometry catalogs, which are combined to cre- 
ate catalogs of suhstellar candidates and used to identify 
targets for extensive spectroscopic follow-up. Our obser- 
vations are designed to reach mass limits of 0.005 Mq 
and below, with the specific goal of probing the bottom 
of the IMF. So far, we have published results for three 
regio ns: NGC 1333 (|Scholz et al.||2009l [20 1l|, p Ophi 
uchi ( Geers et al.||20rT ), and Chaniaeleoii-1 ( Muzic et al 
Mil). 



p Oph is one of the closest regions with active star 
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formation. Its distan ce is estimated to be (125 ± 25) pc 
( [deGeuset al.|1989[ ). The main cfoud, LI 688, is a dense 
molecular core" with visual extinction up to 50 ~ 100 mag 
(Wilking & Lada 1983). It hosts an embedded infrared 
cluster ot around 200 stars, inferred to have a median 
age of 0.3 Myr, and surrounded by muhi ple clusters of 
young stars with a median age of 2.1 Myr (Wilking et al. 
2005 and references therein). Currently, there are ~ 25 



spectroscopically confirmed BDs in p Oph (Wilki ng et al. 
^0081 lAlves de Oliveira et al.||2010| [MaFsh et al 



Geers et al.|201ip . The current census is, however, likely 



incomplete due to strong and variable extinction. When 
compared to the other two SONYC targets NGC1333 
and Cha-I, p Oph suffers much stronger extinction (a 
few tens of magnitudes vs. Ay < 10 mag), p Oph is 
not as compact as NGC1333, but is more compact than 
Cha-I. Finally, Cha-I is slightly older than the other two 
clusters (2-3 Myrs vs. < IMyr). 



In our earlier paper on p Oph (Geers et al. 2011) we 
spectroscopically confirmed a previously unknown brown 
dwarf, and identified a number of new candidates based 
on Spitzer data. In this paper we present the results of a 
second round of spectroscopy, based on our photometric 
catalogs and candidate selection presented in the Geers 
et al. paper. We have obtained about 100 new spectra 
in L1688, using the multi-object spectrograph FMOS at 
the Subaru Telescope, and the integral-field spectrograph 
SINFONI at the ESO's Very Large Telescope (VLT) 
This paper is organized as follows. Photometry and cri 
teria for candidate selection are explained in Section 
and the spectroscopic follow-up in Section [3] The resu 
are presented in Section [4] and discussed in Section 
Conclusions are summarized in Section [6l 

2. MULTI-BAND PHOTOMETRY AND SELECTION OF 
CANDIDATES 

Observations, data reduction, as well as the cross- 
correlation of the optical, near -infrared and Spitz er cat- 



ts 



alogs, are described in detail in [Geers et al. (2011 ). Here 
we give a summary of the facilities and filters used for 
the observations, and the description of the photomet- 
ric catalogs used to select candidates for the follow-up 
spectroscopy. 

The optical images in the SDSS i' filter were obtained 
with the Subaru P rime Focus Camera ( Suprime-Cam) 
wide field imager ( [Miyazaki et al.||2002[ ). The Multi- 
Object Infrared Ca mera and Spectrograph (MOIRCS; 
Suzuki et al.| 2008) at the Subaru telescope was used 
to observe p Oph m the J- and Ks-hands. To increase 
the dynamic range of the survey in the near-infrared, we 
retrieved additional JKs photometry from the 2MASS 

We also re- 



point-source catalog ( Skrutskie et al 



2006) 



trieved the Spitzer IRAC and MIPS photometry of p Oph 
from the Legacy Program data archive available at the 
Spitzer Science Center, using the High reliability cata- 
log (HREL) created by the Cores to Disks (c2d) Legacy 
team. The 3.5tT detection limit of our survey is at i' = 26, 
J — 22, and Ks — 19.5 mag. The completeness lim- 
its correspond to i' = 24.15 ± 0.3, J = 20.6 ± 0.3, and 
Ks = 17.8 ± 0.3 mag. 

To identify candidates for spectroscopy, we use two dif- 
ferent selection procedures: 

(1) We start with the catalog of common sources ap- 
pearing in i, J, and Ks catalogs. Low-mass and BD 



candidates are selected as objects detected in all three 
bands, and exhibiting a color excess in the {i' , i' — J) 
color-magnitude diagram. These are objects located to 
the right of the candidate selection cutoff line in Figure[2j 
This iJKs selection results in a catalog of 537 candidate 
sources shown as dots in Figure [T] 309 of these sources 
come from the MOIRCS JKs catalog (black circles in 
Figure [2]), and the rest from 2MASS (purple circles); and 
(2) We cross-match JKs catalogs with the Spitzer cata- 
log and select objects with IRAC colors consistent with or 
redder than sources with a circumstellar disk (Figure ro| . 
J K / Spitzer selection yields 83 candidate members of p 
Oph, within the boundaries of our optical survey. These 
objects are marked with '-f' in Figure |T] 26 of these are 
found to be BD candidates based on their near-infrared 
magnitude and color, with 12/26 designated as spectro- 
scopically confirmed BDs in the literature. 

35 out of 83 candidates from the JK/ Spitzer list have 
an optical counterpart. 33/35 are also identified as can- 
didates in the iJKs selection, while the remaining two 
are located to the left of the candidate selection line in 
Figure [2] We obtained a spectrum for one of these two 
sources, and found that it is probably not a member of 
p Oph. We note that the majority (32) of the 35 over- 
lapping JK/ Spitzer candidates are part of the 2MASS- 
based iJKs catalog, and only 3 are from the MOIRCS 
iJKs candidate selection. 

3. FOLLOW-UP SPECTROSCOPY 

In our first spectroscopy run in May 2009, we used 
MOIRCS/Subaru to obtain the spectra of 58 BD candi- 
dates selected from our iJKs c atalog, and one p reviously 
known BD candidate member ( Geers et al.|2011|) . The 58 
selected candidates all come from the iJ Ks catalog based 
on the MOIRCS JKs photometry (309 objects). Dur- 
ing 2010 we were awarded time at the Subaru Telescope 
and the VLT that allowed us to obtain spectra of ~100 
additional substellar candidates in p Oph. These new ob- 
servations, data reduction and number statistics of the 
observed candidates from the two catalogs are presented 
in this section. Candidates with spectroscopic follow-up 
are marked with crosses in Figure [2] and squares in Fig- 
ure El 

3.1. FMOS/Subaru 

We obtained spectra in the near-infrared J- and H- 
bands (0.9 - 1.8 pm) using the Fiber Multi Objects 



Spectrograp h (FMOS) at the Subaru Telescope (iKimura 

[US 



et al. 2010). The fiber positioner system of FMub con 
figures 400 fibers in the 30-arcmin diameter field of view, 
and the spectra are extracted by the two spectrographs 
(IRSl and IRS2). Our data were obtained during the 
SlOA semester in a shared-risk mode, using only one of 
the spectrographs (IRSl; 200 fibers) in the low-resolution 
mode (R ~ 600). Data reduction was carried out using 
the software package supplied by Subaru during the ob- 
serving run. The package consists of IRAF tasks and C 
programs using the CFITSIO library. The recommended 
telluric standards for FMOS are stars of spectral type F, 
G or early K, at a known (and preferably low) extinction. 
As the stars in our science field suffer strong extinction, 
we observed an additional field at the same average air- 
mass as p Oph, containing several GOV stars that were 
used for telluric correction. 
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Fig. 1. — S patial distribution of sources in p Opli. Contours are Av = 5, 10, 15, 20, 25, 30, as derived from 2MASS by the COMPLETE 
project I Ridge et al.| |2006: available at www.cfa.harvard.edu/COMPLETE). Dotted line: j'-band imaging coverage; dashed line: J- and 
i^5-band imagmg coverage; the circle shows the FMOS field- qf-view. Carididat e members are indicated with dots {i'JK selection), and 
pluses (J K I Syitzer selection). SON YC spectroscopic targets | |Geers et al.|20lT| and this work) are marked with crosses, p Oph members 
from[Wilking ct al . (2008) and Eric kson et al.| ( |2011| l are shown as blacic diamonds, while the known BDs | |Wilking et aL|[2008| |Alves de| 



|01iveira~et al.|2UlU| jMarsh et al.|2Dlljj | are shown as red diamonds. Red squares mark the BDs and VLM stars tound m the framework ot 
SOM VC. 



We observed p Oph during the two nights of June 3 — 4 
2010, and obtained in total 81 spectra. 14 of the objects 
came from the J K / Spitzer candidate hst, 46 from the 
iJKs-haseA candidate selection, and the remaining fibers 
were placed on a number of known members of spectral 
type K and M, to serve as a reference. We note that the 
iJKs candidate selection is based only on the MOIRCS 
JKs catalog (309 objects). 

3.2. SINFONI/VLT 

The spectra of 13 objects from the J K / Spitzer can- 
didate list were obtained using the integral-field spec- 
trograph SINFONI at the ESO/VLT, in the H- and K- 
bands (program ID 385.C-0450). Two of the objects have 
also been observed with FMOS. The observations were 
carried out one object at a time, using the largest avail- 
able FOV offered by SINFONI (8" x 8") and in the no- 
AO mode. The spectral resolution delivered by the H-l-K 
setting of the instrument is R ^ 1500. Data reduction 
was performed using the SINFONI pipeline supplied at 
ESO. Telluric correction was performed with the help 
of standard stars observed before and after our targets. 
Standard stars were all of the early B-type, which ex- 
hibit a large number of hydrogen lines throughout the 
H- and iiT-bands. Some of the lines can be manually re- 



moved from the standard-star spectra (e.g. Br 7, Br 10, 
Br 11), while others appear blended with the telluric fea- 
tures and thus cannot be easily removed. Some residuals 
due to uncorrected hydrogen lines appear in our spectra, 
but do not affect the results, as we are interested only 
in broad features. These residuals mainly mimic emis- 
sion lines in the 1.55 - 1.65 pLVn region in some of the 
SINFONI spectra. 

4. SPECTRAL ANALYSIS 
4.1. Spectral fitting procedure 
The spectral analysis is identi cal to the one p resented 



Scholz et al I ([2009) and ,Geers et al.| ( |2011 ). Spec 
tra ot young objects later than M5 show a cha racter 

p005) 
h sides 



istic broad peak in the if- band ( Gushing et al. 
The feature is caused by water absorption on bot' 



of the iJ-band and its depth depends strongly on ef- 
fective temperature. While the i7-band peak appears 
round in old field dwarfs, it is unambiguously triangular 
in yo ung, low-gravity so urces (e.g. [Brandeker et al.^2006[ 
jKirkpatrick et al.|2006 ). In addition, young BDs are ex- 
pected to have fiat or increasing ifg-band spectra with 
CO absorption bands at A > 2.3^m, and a sharp drop in 
flux at the edge of the J-band (^ 1.3/.tm), again due to 
water absorption. 
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Fig. 2. — («',«' — J) color-magnitude diagram. Filled circles rep- 
resent all sources from our i' JK catalog (MOIRCS JK), with dark 
symbols indicating the BD candidates. Sources from the 2MASS- 
based i' JK catalog are indicated with purple open circles. Crosses 
show SONYC spectroscopy follow-up targets, with spectroscopi- 
cally confirmed VLMOs marked with squares. Among SONYC 
VLMOs, the object found at («' - J,i') = (5.2,24.1) is identified 
as a low-mass star, while the remaining five are BDs. Spectroscop- 
ically confirmed brown dwarfs from the literature fWilki ng et al.| 
^08 Alves dc Oliveira et al. 2010) are shown as diamon ds, while 
stars m dicate higher mass members of p Oph from Wilking et al.' 
l|2008[ l andlErickson et al. (2011). Atmosphere model isochrones 
tor age 1 Myr are ovcrplottcd tor DUSTYOO ( |Chabrier et al..2000j 
and COND03 (Baraffc et al. 2003). The solid lines denotes the tiU 
candidate selection cutott. The g ; reen dash-dotted line denotes the 
new selection cutoff (see Section]5.2|for explanation). 



The spectra were first assessed through visual inspec- 
tion. Very noisy spectra 5% of all spectra) and those 
lacking the above mentioned features characteristic for 
BDs (~ 75% of all spectra) are discarded. The latter 
sources exhibit smooth, featureless spectra, which make 
it difficult to constrain their nature. The majority of 
these objects are expected to be background late-type 
giant or dwarf stars. Some of them could be embedded 
YSOs with spectral types earlier than M, although we do 
not have any evidence for youth. This selection yields 21 
VLM candidates. 

In a second step, we estimate the effective tempera- 
tures of the selected VLM candidates by fitting model 
spectra from the DUSTY series (Allard et al. 200 1J, fol - 
lowing the procedure outlined in jSchoiz et al.| ( j2009[ ) . 
Since extinction and effective temperature are difficult 
to determine separately with low-resolution, low signal- 
to-noise spectra, we started by adopting the Ay as de- 
termined from the J — Ks colors, assuming an intrin- 
sic BD photospheric J — Ke; color of 1 mag. We apply 
the extinction law from Cardelli et al. (1989), assuming 



Rv — 4. We vary the effective temperature from 1800 K 
to 3900 K, in steps of 100 K, keeping the extinction at a 
fixed value. The best fit temperature is determined from 
a analysis. Although in almost all cases this procedure 
results in a decent match between the observed spectra 
and the best-fit models, we note that small adjustments 
to the value of the extinction (AAv < 2) can slightly 
improve the quality of the fit. This adjustment is de- 
termined by changing Ay in steps of 1 around the value 



derived from the J — Ks color, rounded to the closest 
integer. For 21 out of 23 entries in Table [l] either a zero 
or a small adjustment (up to Av=2) is required, while 
in the remaining two cases AAy 3 results in a better 
fit. Both values of Ay are listed in Table [T] 

As demonstrated in Scholz et al. (2009j), it is often dif- 
ficult to reliably distinguish between two models with 
AT = lOOi^:. However, for AT > 2mK there are 
discernible differences between models and observations, 
combined with a clear increase in x^. Thus, a conserva- 
tive estimate for a typical uncertainty in Tcff is ±200 K. 
In cases where an equally good fit is obtained using two 
models separated by 100 K, as the best-fit T^i in Table [T] 
we list the average value, while in Figures^ and |4] we plot 
the model with the lower T^s- Although the temperature 
fitting grid extends up to 3900 K, we note that our fitting 
scheme does not work that well above 3200 K. The rea- 
son is the shape of the H-band peak that becomes fiat, 
making it difficult to discriminate between the models 
at higher temperatures. The situation is further compli- 
cated by high extinctions typically encountered in the p 
Oph region. We find spectra of three objects to be consis- 
tent with effective temperatures of 3200 K or somewhat 
above. We set 3200 K as a lower temperature limit for 
these objects in Table [T] 

4.2. Results of the spectral fitting 

The spectra of 21 BDs and very-low-mass objects 
(VLMO) are shown in Figure |3] (FMOS) and Figure jl 
(SINFONI). Two objects appear in both datasets (#T 
and #11), The best fit Toff resulting from the two inde- 
pendent fits are consistent, differing by only 50 K and 100 
K, respectively. Spectra are shown in black and the best- 
fit models in color. For clarity, the models are slightly 
offset from the spectra. Both models and spectra have 
been smoothed with the boxcar function of 30 A prior to 
plotting, to make it easier to appreciate the quality of the 
broad-band fit. The 21 BDs and VLMOs with SINFONI 
and/or FMOS spectra are listed in Table [T] including the 
best model fit parameters. 

Our observations result in 13 probable substellar mem- 
bers (defined here as objects with T^ff < 3000K) and 8 
very-low-mass stars (3000 K < Teg < 3200K). Among 
these objects, we report 6 new likely BDs and 2 new 
very-low-mass stars. We label those objects SO NYC- 
RhoOph-3 to 10, following the nomenclature of Geers 
et al. ( |2011 ) who published the first two discoveries. The 
SIX new BDs have effective temperatures between 2600 K 
and 3000 K, which, according to the IMyr DUSTY and 
COND isochrones, correspond to the masses between 
0.03 Mq and 0.1 Mq. 

The H-band spectrum of the object #7 (SONYC- 
RhoOph-5) appears to have a slightly rounded shape 
when compared to the spectra of other objects in Fig- 
ure [31 However, neither the comparison with the spectra 
of field dwarfs or giants provides a satisfactory match 
with the observed spectrum (see Figure 5]). The source 
is detected in all four IRAC bands anci shows a clear 
mid-IR excess: with 11-12 = 0.41 and 13-14 = 1.03, it 
resides within the class-II object area in Figure |6] We 
thus conclude that SONYC-RhoOph-5 is a probable new 
BD member of the cluster. 

The H-band peak around 1.68/im in the spectrum of 
the object #14 might not appear as sharp as expected 
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Fig. 3. — Spectra of yery-low-mass obje cts in p Oph, observed with FMOS/Subaru. The spectra are shown in black, with the best-fit 
models (AMES-Dusty; [AUard et al.|2001^ in color. The models are slightly offset for clarity. The spectra and the models were smoothed 
prior to plotting, so that the quahty ot the fit to the broad-band features can be appreciated. The models have been reddened by the 
best-fit Ay, given in Table [l] The spectra are plotted in the order of increasing Teff from bottom to top and from left to right. 
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Fig. 4. — Spectra of very-low-mass objects in p Oph, observed with SINFONl/VLT. See caption to Figurejsjfor plot details. 
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Fig. 5.— Spectrum of the object #7 (black; SONYC- 
RhoOph-5) in comparison with the best-fit low-gravity model 
(blue) and spectra of three Jield dwarfs of spectral type M5 
(2MASSJ015 32750+3631482: Burgasscr et al. 2004), M7 ( VB 8; 
Burgasser et al. 2008 ) and M9 (2MASS J 12425U52+2357231 |Kirk | 
Patrick ct al. 201U), as well as a spectrum of an M8 giant (tiU 
TT3285; bpcA Libraries) shown in grey. The model and template 
spectra have been reddened by Ay = 14 and normalized at 1.65/xm. 



for a young BD with Tcs — 2550K, but this can be at- 
tributed to the strong no ise in this spectrum. In |Alv es| 



de Oliveira et al. (2010) one can clearly see that this 
source exhibits the sharp H-band peak characteristic for 
young BDs. They classify the source as M8.25, which is 
consistent both with the effective temperature obtained 
from our fits, and t he sp ectral type derived from various 



indices (see Section 4.3| and Table [2|). 

In this paper and in previous papers of the SONYC 
series, we have chosen to draw a line between the stel- 
lar and substellar regime either at T^s = 3000 K or at 
SpT = M6. According to the COND03 and DUSTY mod- 
els, a 1 Myr old object at the substellar boundary will 
have an effective temperattire of about 3000 K. In the 



scale defined by Luhman et al 
responds to a spectral type M6 



( |2003[ ), this limit cor- 
However, the relation 



between the effective temperature and spectral type for 
young VLMOs is not that well established. There are two 
objects that have at least once been classified as brown 
dwarfs in the literature (M6 or later) and that we do 
not define as substellar. The first o ne is the source #21 
whi ch we find to have Teg > 3 200K. |Wilking et al] ( |1999 1 
and Luhman & Rieke ( 1999 1 classify this source as M7 
and M6.5, respectively. Based on our fits, we identify 
this source a very- low-mass star, rather than a BD. 

The second object where we find discrepancy 
with the literature classification is found at (a, 
^)20oo =(16:26:31.36, -24:25 :30.2) and was classified as 
M6 by [Wilking et al.|(|1999|) from th e NIR data, but also 
as M5 (INIR; [Luhman fc Rieke|1999p and M3.75 (optical, 



much tiatter than in BD spectra, resulting in T^s clearly 
above 3200 K. Thus, we discard this object as a potential 
BD, and conclude that our spectrum is more consistent 
with the earlier spectral classification from the literature. 

4.3. Spectral type indices 

In this section we attempt to assign spectral types to 
the VLMOs from our survey, based on various indices 
suggested in the literature. The SINFONI spectra are 
particularly suited for a comparison of different methods, 
thanks to the spectral coverage, resolution, and overall 
quality of the spectra. We determine spectral types based 
on the following three indices 



( [2011 1 
mixes m 



(1) H-peak index (HPI) defined by Scholz et al 
and based on the ratio between the average 
1.675 — 1.685/im and 1.495 — 1.505^m wavelength re- 
gions. The HPI takes advantage of the slope on the 
blue side of the H-band peak observed in late-M and 
L-type objects, and is related to the spectral type as 
SpT= -0.84 + 7.66 -HPI. The index is highly sensitive 
to the spectral type in the >M7 regime, but cannot be 
used reliably for earlier spectral types where the slope of 
the feature becomes too shallow. 

(2) Q-index defined in Wilking et al.[ ( [1999[ ) as Q = 
{Fl/F2){F3/F2y. Fl, F2, and F3 are the average val- 
ues of the relative flux densities computed in the three 
bands 2.07 - 2.13, 2.267 - 2.285, and 2.40 - 2.50/im, re- 
spectively. The useful portion of our SINFONI spectra 
extends only to ~ 2.45/im. As evident from Figure [4j the 
continuum of the spectra longward of ^ 2.3/im can be ap- 
proximated by a linear function. We fit a linear function 
to the wavelength range 2.30 — 2.45/im (2.43/im in case 
of the object #4) and extrapolate to the 2.45 — 2.50/im 
region. We use (3 = 1.26 and the SpT — Q relation de- 



rived for p Oph by Gushing et al. ( 2000 
(3) The iJsO-index proposed by [Allers et al.| ([2007 



Eg. 1 and 2) 



Eq. 1). This index is defined as the ratio between the 
average fluxes in 1.550 — 1.560/im and 1.492 — 1.502fim 
wavelength regions, and is linearly dependent on spectral 
type between M5 and LO. 

All spectra have been extinction corrected using the best- 
fit value of extinction obtained from our model fits (see 
SectionEl), a nd assuming the extinction law from Cardelli 
et al. (lS89|. The spectral types resulting from the three 



indices are shown in Table ^ The wavelength coverage 
of the FMOS data does not extend to the K-band, and 
thus they cannot be used to calculate the Q-index. For 
the object #14, one of the narrow band regions used to 
calculate the lf20-index coincides with a strong noise 
feature, which prevent us to use this index to derive a 
reliable spectral type. For the objects #4 and #11, we 
list spectral type estimates based on their FMOS and 
SINFONI spectra, in that order. 

We can make a comparison between the spectral types 
derived from different indices, excluding those where we 
only list the limits on the spectral type. The quoted 
uncertainties for the spectral typ es derived f rom the three 
indices are ±0.4 for HPI (Sch olz et al.|2011[ ), ±1.0 for the 
HaO-index (Allers ct al. 2007j), and ±1.5 tor the Q-index 
( [Wilking et al.||1999) . When comparing the HPI-derived 
values with those based on the other two indices, we find 
that our spectral type scheme gives systematically later 
spectral types (in all cases except one). On average, the 
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TABLE 1 

Parameters of p Oph M-type objects observed with FMOS/Subaru and SINFONI/VLT 



4k 

Tr 


a(J2000) 


(5(J2000) 








CIT V / 


SpT 






1 


16 26 03.28 


-24 30 25.8 


F 


8.4 


10 


2650 






SONYC-RhoOph-7 


2 


16 27 05.93 


-24 18 40.2 


F 


7.8 


8 


2700 






SONYC-RhoOph-6 


3 


16 27 38.63 


-24 38 39.2 


F 


6.6 


8 


2750 


M8.5, M7, M6 


123 


GY310 


4 


16 26 22.27 


-24 24 07.1 


F 


9.7 


8 


2750 


M6.5, M8.5 


13 


GYll 


4 


16 26 22.27 


-24 24 07.1 


S 


9.7 


12 


2800 


M6.5, M8.5 


13 


GYll 


5 


16 26 18.82 


-24 26 10.5 


F 


9.4 


10 


2850 


M7.5, M5.5, M7 


123 


CRBR 14; ISO-Oph-23 


6 


16 26 18.98 


-24 24 14.3 


F 


17.6 


16 


2850 


M5 


12 


CRBR 15 


7 


16 26 22.96 


-24 28 46.1 


F 


13.2 


14 


2900 






SONYC-RhoOph-5 


8 


16 26 51.28 


-24 32 42.0 


F 


2.3 


2 


2900 


M8, M8.5 


54 


GY141 


9 


16 26 53.35 


-24 40 02.3 


F 


12.4 


13 


3000 






SONYC-RlioOph-3 


10 


16 26 22.19 


-24 23 52.4 


F 


13.8 


14 


3100 


M8.5, M6.5 


12 


GYIO 


11 


16 26 27.81 


-24 26 41.8 


F 


6.4 


6 


>3200 


M5, M6 


61 


GY37 


11 


16 26 27.81 


-24 26 41.8 


S 


6.4 


6 


3100 


M5, M6 


61 


GY37 


12 


16 26 18.14 


-24 30 33.0 


F 


14.7 


13 


3200 






SONYC-RhoOph-4 


13 


16 27 46.29 


-24 31 41.2 


F 


8.3 


9 


>3200 


M6 


3 


GY350 


14 


16 27 40.84 


-24 29 00.7 


S 


3.1 


3 


2550 


M8.25 


7 


CFHTWIR-Oph96 


15 


16 26 18.58 


-24 29 51.4 


S 


15.0 


18 


2600 






SONYC-RhoOph-8, CFHTWIR-Ophl6 


16 


16 27 26.58 


-24 25 54.4 


s 


0.9 


2 


2650 


M8 


6 


GY264 


17 


16 27 26.22 


-24 19 23.0 


s 


13.6 


15 


2800 






SONYC-RhoOph-10: CFHTWIR-Oph78 


18 


16 26 37.81 


-24 39 03.2 


s 


6.2 


6 


3100 






SONYC-RlioOph-9; CFHTWIR-Oph31 


19 


16 27 06.60 


-24 41 48.8 


s 


3.7 


6 


3150 


M5.5, M6 


63 


GY204 


20 


16 26 23.81 


-24 18 29.0 


s 


8.2 


10 


3200 


M6.7, M7 


58 


CRBR 31 


21 


16 27 05.98 


-24 28 36.3 


F 


14.7 


17 


>3200 


M7, M6.5 


12 


GY202 



Note. — Based on AMES-Dusty models jAllard et al.|2001[ 
Instrument used for spectroscopy. F: FMOS/Subaru, S: SINFONI/VLT 



Visual extinction in magnitudes, and calculated from J 
Best-fit extinction from spe ctra.! fittinjy^ in yna.jynj tuHpj^ 
Spcctral-typc references: [1] 



K color, assuming the intrinsic J — K oi\. 



trfiil ntting, in.CTiaigii)tud(^ 

Wilking et al.|l|l999|l; [2]|Luhman &: RiekeHl999^; [3]|Natt a_et al.H2002[ ; [4] |Luhman et al.|jl997t ; [5] |Cushing et al.| 



|Bontemps et al.| | [2001^ . 



2000[ ; [6] |Wilking et al.| l [2005t ; | 7||Alves de Oliveira et al.M2010| |8| IWiiking 
References for identi Sers: GY: Greene fc Young| ]1992[ ; CFHTWIRT l^lves 



et al. ( 200 8jr 

cle Olivcira et al.| l [2oTo) ; CRBR: [Comeron et al] l |1993[ l; ISO-Oph: 



TABLE 2 

Spectral types calculated from the SINFONI 
spectra 



# 


HPP 


Spectral Type 


H20'^ 


1 


M9.7 




M9.6 


2 


M8.6 




M6.9 


3 


M8.3 




M6.6 


4 


M8.5, M8.0 


M6.1 


M6.6, M7.8 


5 


M8.5 




M8.0 


6 


M8.0 




M5.7 


7 


M7.4 




M5.7 


8 


M8.0 




M6.3 


9 


M7.5 




M5.1 


10 


M7.4 




M5.5 


11 


<M7, M7.1 


M5.8 


<M5, M5.8 


12 


M7.0 




M5.0 


13 


M7.1 




<M5 


14 


M9.3 


M7.1 




15 


M9.1 


M8.7 


M8.0 


16 


M8.6 


M8.0 


M8.0 


17 


M7.8 


M6.2 


M8.4 


18 


M7.0 


M6.8 


M5.4 


19 


M7.1 


M5.9 


M6.4 


20 


M7.0 


M6.6 


M5.7 


21 


<M7 




<M5 



^ Derived from the HPI ind ex jScholz et al.|201l|| 

Der ived from the Q-index jWilking et al.|l999]rCushing et al.l 

12000} . 

Derived from the HaO-index l |Allers et al.|2007[ 

difference between the HPf and H2O spectral types is 1.3. 
The difference is more pronounced for the FMOS data: 
the average value is 1.6, while for the values derived from 
the SINFONI data we get better agreement (0.9 spectral 
types). This can probably be attributed to the better 



spectral resolution and signal-to-noise of the SINFONI 
data. The average difference between the HPI and Q- 
index spectral types is 1.1. 

We tested the effect that uncertainties in the fitted 
Ay have on the derived spectral types. A change in 
the extinction of A^y = 2 results in AS'pT ~ 0.6 sub- 
types, for the spectral types derived from the HPI and 
the iJ20-index, while the Q-index is esentially reddening- 
independent. The change in the Ay affects the HPI and 
HiO spectral types for our sample objects in a similar 
way. A possible explanation for the discrepancies be- 
tween spectral types derived from the HPI and other 
two indices might be in the extinctions of the sample 
from which HPI was derived. HPI was derived for ob- 
jects with low extinction, i.e. the extinction correction 
was less problematic. The majority of our objects in p 
Oph are at large extinctions, which will introduce ad- 
ditional uncertainty to the spectral type as derived from 
the HPI. On the other hand, the sample that was used to 
derive the Ii20 index covers a larger range of extinctions. 

Given the uncertainties, we conclude that the spectral 
types derived from the different indices are in agreement, 
with the HPI spectral types being systematically later by 
about one spectral subtype from the values derived from 
the other two indices. 

5. DISCUSSION 

5.1. Objects with mid-infrared color excess 

Figure [6] shows a color-color diagram constructed from 
Spitzer IRAC photometry. The dashed line denot es the 
area where class II objects are located, based on Allen 
et all (|2004|). For 83 objects, the IRAC colors are con- 
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Fig. 6.— Left: Color - color diagram in ([3.6] - [4.5], [5.8] - [8.0]), constructed from Spitzer IRAC photometry, for 2 MASS and MOIRCS 
sources with a Spitzer match (black dots). The dashed line denotes the area where class II objects are located, based on jAUen et al.||2004[ l. 
Spectroscopically confirmed members of p Oph from Wilking et al. (2008) and Erickson et al. (2011) are represented by diamonds, with the 
confirmed brown dwarfs (defined as the objects with spectral type > M6, from Wilkmg e t al.|20U8| |Alves de Oliveira et al.|2010l and this 
work) shown in red. Right: Same as left; open circles mark the objects with spectra obtained as part of SUIN Y(J (bve usmg MUIkOS/MOS, 
and the rest using FMOS and SINFONI). SONYC BDs and VLMOs are shown as the red symbols. 
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have been either observed by us, or have been spectro- 
scopically confirmed as members in the literature. 15 of 
them have spectral type later than M5 or estimated effec- 
tive temperature at or below 3000 K. Thus, we conclude 
that the number of missing brown dwarfs in this sample 
is (83 — 65) X 15/65 w 4. Of course, the total number 
of missing substellar objects in p Oph is probably much 
higher. Our estimate is valid only for the objects with 
disks (i.e. exhibiting MIR excess), and down to a certain 
limit in mass and extinction. This limit is primarily im- 
posed by the sensitivity of the Spitzer data. The solid 
line in Figure [7] shows the approximate limit of the com- 
bined J K / Spitzer catalog. It corresponds to ~ O.OIM© 
at A y = Oand~ O. OSMfD at Ay = 30, based on COND03 



4 
J-K 



Fig. 7. — Near-infrared color-magnitude diagram for the sources 
in p Oph, based on MOIRCS and 2MASS photometry. Spectro- 
sco pically confir med stellar and substellar me mbers ( Wilki ng et al.| 
[m)8, .Erickson eFal^lMTT] [Aives de Oliveira et al. ,2010, and this 
work) are marked with crosses; objects with Spitzer counterpart 
are shown in red. The solid line shows the approximate limit of 
the combined JK — Spitzer catalog. 

sistent with or redder than sources with a circumstel- 
lar disk, here defined as (13 - 14) > 0.4 and (II - 12) 
> 0.0. 46 of these objects appear in the list of spectro- 



( [Baraffe et al.|2003| and DUSTYOO ( [Chabrier et al.|2000 l 
evolutionary models. 

Based on the ISOCA M survey of the main p Oph core, 
Bontemps et al. (2001) find that the fraction of objects 
with disks among the embedded cluster YSO population 



is ^ 54%. Fr om a small sample of 9 BDs, Jayawardhana 



scopical ly confirmed p Oph members (diamonds; Wilking 
et al.|20 08[). 33 objects have been observed in the frame 



work of'SON YC (black squares). The two samples have 
19 objects in common. Ther e are 12 spectroscopi cally 
confirmed brown dwarfs from [Wilking et aL ( 2008 ) and 



Alves de Oliveira et al. (20101) and 3 identified in this 



work appearing in this plot. I'here are 18 objects with 
Spitzer excess that, to our knowledge, do not have avail- 
able spectra (objects shown only as crosses to the right 
of (13 - 14) 0.4 in Figure [6|. 

Of 83 candidate members with mid-infrared excess, 65 



et al.| ( 2003 ) calculate the disk fraction among the BD 
population of ~ (67 ± 27)%. The number of the missing 
BDs with disks in our sample is ~ 4. Accounting for 
the disk fraction, we estimate that the total number of 
missing BDs ^ 7. 

5.2. Substellar population in p Oph 

Based on the results of our survey, and combined with 
literature data, we can try to put limits on the total num- 
ber of substellar objects in p Oph. However, we must 
note that p Oph, with its high and strongly varying ex- 
tinction, is a difficult field to search for BDs. All current 
surveys, including SONYC, have to be considered incom- 
plete. 

SONYC was designed to mitigate selection biases that 
are typically induced by commonly used mid-infrared- 
based surveys (sources with discs), or narrow-band 
methane surveys (T-dwarfs). Our primary means of se- 
lecting candidates for spectroscopy is broad-band optical 
and near-infrared imaging that aims to detect the photo- 
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sphere. This selection is clearly not effective in the areas 
with high extinction, and is also significantly affected by 
the background contamination. However, in an attempt 
to get an estimate of the number of missing substellar 
objects in our survey area, without being biased towards 
the disk objects, we base the following discussion only on 
the optical-NIR candidate selection. 

Figure [2] shows the (i', i' — J) color-magnitude dia- 
gram for all the sources from our i'JK catalogs. Black 
circles show 309 candidates from the MOIRCS JK se- 
lection, while purple open circles show sources with 
JK photometry from 2MASS. Crosses show SONYC 
spectroscopic targets, with confirmed VLMOs shown as 
squares. Among the SONYC VLMOs, the object found 
at — J) = (24.1, 5.2) is identified as a low-mass star, 
while the remaining five are BDs. We also show spec- 
troscopically confirmed stellar members (stars) and BDs 
(diamonds) from the literature. Distribution of the ob- 
jects with spectra and confirmed substellar member can- 
didates shows that the criterion used for the candidate 
selection is fairly conservative. 

Candidates from the optical-NIR selection that have 
been spectroscopically followed up come from the i' JKg 
catalog based on the MOIRCS JK photometry. This 
candidate list contains 309 objects. We have obtained 
spectra for 102 of these objects. Five were confirmed as 
likely BDs. Thus, the spectroscopic follow-up suggests 
that our catalog could contain (309 - 102) x 5/102 sa 10 
additional BDs. This number might actually be smaller 
because, as evident from Figure^ the confirmed BDs 
tend to be located to the right ofthe DUSTYOO 1 Myr 
isochrone. Based on this observation it might be reason- 
able to shift the selection cutoff line slightly towards the 
right (green dash-dotted line) . The probability of finding 
any new members in the strip between our old selection 
line and the new one is expected to be very low. The 
new selection cutoff was chosen to lie approximately at 
the same distance from the two isochrones shown in Fig- 
ure [2j The new candidate list then contains 151 objects, 
75 of which have spectra. The expected number of miss- 
ing BDs is then (151 - 75) x 5/75 w 5. 

The combined i' JK catalog is shallower than the i'- 
band catalog and its completeness limit is at i' « 23.4. 
The majority of the optically selected candidates lie in 
the area with low extinction (see Figure [T 1 : 89% of the 
candidates are located at Ay < 15, as judged from the 
COMPLETE map . However, as shown in Figure 5 of 
Geers et al.| ( poTlj ), all objects from the MOIRCS i'JK 
catalog, with the exception of only two sources, prob- 
ably suffer optical extinctions of less than 15. The ex- 
tinction for each source was calculated based on the ob- 
served (J — K) colors, assuming an intrinsic {J — K) of 1. 
Thus, even in the regions with apparently large extinc- 
tion, our optical-NIR observations seem to be penetrat- 
ing the cloud only down to a certain depth, corresponding 
to Ay of about 15. The completeness limit of i' « 23.4 
corresponds to masses 0.003 — O.O3M0 for Ay =0—15, 
according to the DUSTYOO and COND03 models for 1 
Myr old objects at the distance of p Oph. We conclude 
that the number of missing objects in our survey area 
down to this mass and extinction limit is 5-10. This num- 
ber should be added to the number of missing o bjec ts es- 
timated from the JK — Spitzer data in Section 5.1 The 
analysis presented in this section is based on the 309- 



candidate MOIRCS i'JK catalog and does not include 
the 2MASS i'JK selection. The JK- Spitzer candidate 
selection is, on the other hand, almost entirely based on 
the 2MASS JK photometry (72/83 candidates). In Fig- 
ure [7| the JK — Spitzer candidates are mostly found 
above the grey fine, while ~ 85% of the MOIRCS i'JK 
candidates lie below it. The two catalogs overlap for only 
3 objects. 

In summary, we conclude that the number of missing 
substellar objects in our survey area is 12-17. The actual 
number of missing substellar objects is probably higher, 
as we certainly might be missing objects at lower masses 
and higher extinctions. 

5.3. Ratio of stars to brown dwarfs 

Based on the currently available data we can put con- 
straints on the ratio of the number of l ow-mass stars to 
the nu mber of BDs in our survey area. [Andersen et al.| 
(20081 use a mass range m — (0.08 — 1}Mq for stars 
^Hd^ = (0.03 - O.O8)M0 for BDs; we call this ratio 
Other authors use m — (0.08 — 1O)M0 for stars 
and m = (0.02 - O.O8)M0 for BDs, hereafter called R2 



dBriceiio et al.||2002l IMuench et al.||2002| iLuhrnan et al 

^' " ^' " 

The two most complete estim ates of the stellar popu- 
lati on in p Oph are thos e from Bontemps et al. ( 2001 1 
and iWilking et al | ( [20081 ) ■ |Bonteinpset al.| ( |200 1l) report 
123 Class II YSOs, 95 of which are located within the 
borders of our optical survey. From their Figure 8, we 
infer --80 YSOs in the mass range (0.08- 1)Mq, and ~84 
in the mass range (0.08 — 1O)M0. We scale these num- 
bers with the fraction of all sources located in our FOV 
(95/1 23), and correct for the binary fraction of about 
29% (Ratzka et al. 20051. The number of stars in our 
FOV IS then ~ 80 and ~ 84 for mass ranges (0.08- 1)M0 
a nd (0.08 — 10)Mi^, res pectively. From the member list 



Wilking et al. (20081 we discard objects with spectral 



type M6 or later, restrict the coordinates to match the 
area covered by our survey, add the three VLMOs discov- 
ered in the framework of SONY C, and four stell a r-mas s 
members recently identified by Erickson et al. (2011). 
We count 86 low-mass stars with spectral type estimates 
between G2 and M6 (i.e. with masses between 0.08 and 
1 Mq, and 90 stars in the mass range (0.08 — 10) A/0. 

There are currently 26 spectroscopically confirmed 
probable BDs within the borders of our survey (20 pre- 
viously known and 6 found in this work). We calculated 
absolute J-band magnitudes for the sample of BDs using 
the A y estimated fro m J — K colors, reddening law from 



Cardelli et al. (1989), and by assuming the distance of 
125 pc. We then compare the results with the COND03 
and DUSTYOO 1 Myr isochrones. There are 21 objects 
with m > O.O3M0, and 23 with m > O.O2M0. Among 
these there are 8 objects for which we obt ain mass 
estimates ab ove O.O8M0; 7 come from the list of Wilking 
et al.| ( |2008[ ) and one is our object SONYC-RhoOph-5. 
All of these objects have spectral type around M6, 
i.e. they lie around the boundary separating low-mass 
stars from BDs. Considering that the mass estimates 
are subject to a number of uncertainties (arising from 
the uncertainties in cluster distance, age, extinction, 
photometry, effective temperature, and the models 
themselves), all of these 8 objects could in fact be either 
low-mass stars or BDs. Thus, for the calculation of the 
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ratios Ri and R2, we assume the following numbers: 

(1) stars with m = (0.08 - 1)Mq: 80 - 94, 

(2) stars with m = (0.08 - 1O)M0: 84 - 98, 

(3) BDs with m = (0.03 - O.O8)M0: 13 - 21, 

(4) BDs with m = (0.02 - 0.08)Mo: 15 - 23. 

To calculate the ratios i?i and R2 we adopt the 
mean values of the intervals specified above. A reason- 
able indication of statistical uncertainties associated 
with these numbers is provided by Poisson statistics, 
i.e. we assume the standard deviation a ~ ^fji. where [i 
stands for the mean of each interval. 

Finally, we arrive at i?i = 5.1 ±1.4 and R2 = 4.8±1.2. 
This resul t is in line with t he values from the literature 
(3.3 - 8.5; Andersen et al.|[2008), although higher t han 



the value observed in NGC1333 {R\ = R2 = 2.3; Scholzl 
et al.||2011[ ), and in the ONC (i?i = 2.4; Andersen et aT] 
zOlip . -Note that the values quoted here and in Andersen 
et al. (2008) should be treated as upper limits. The total 
number of BDs is certainly higher than the values given 
above due to the lower mass cutoff set at 0.02 — O.O3M0 
and the missing BDs. 

5.4. Spatial distribution and disks 

In Figure[8]we show the projected distribution of spec- 
troscopically confirmed stellar members and BDs in p 
Oph. Stellar-mass m embers ( black cross es) are from 



Wilking et al.| ( |2008D and [Erickson et al | ( |2011[ ), with 
three additional objects identihcd as very-iow-mass stars 
in the framework of SONYC. R ed crosses show spec - 
troscopically confirmed BDs from 'Wilking et al.^ ('2008| 



Alves de Olive ira et al.|(2010), Marsh eFarT pOlO), Geers] 
et al. I2OII ), Erickson et al.| ( |201 1) and this work. Both 
catalogs of stellar and substellar-mass members were 
cross-matched with the catalog containing objects show- 
ing MIR color-excess, which is interpreted as the evidence 
for an existence of a circumstellar disk around these ob- 
jects. The stellar population clearly shows clustering to- 
wards the inner part of the L1688 core. There are 80 
objects within the boundaries of our optical survey (dot- 
ted rectangle), and 65 objects outside of it (in the area 
roughly twice as large). At a first glance it might seem 
that even stronger clustering exists among the BD pop- 
ulation. However, this can be attributed to the lack of 
a uniform BD survey in the entire cluster region. Most 
of the existing BD surveys, including this one, have been 
concentrated on the inner core of the clou d, and only two 



surveys actually probed t he larger area (Gushing et al. 
Wilking et al.|2005 ). Therefore, the larger star-to- 



2000 



BD ratio in the outskirts of the cluster is probably not 
real. Within our FOV, there is apparently no significant 
difference between the spatial distribution of stellar and 
substellar populations. 

A comparison between the two populations can also 
be made based on the presence of disk excess. Here we 
again concentrate on the area within the i'-band survey 
boundaries. The ratio of the number of stellar mem- 
bers with disks, to the total number of stellar members 
is 37/90 = (41 ± 8)%. This number probably represents 
the lower limit for the disk fraction, because of the rela- 
tively rigorous criteria on photometry that were applied 
to the 2MASS data. For the 2MASS - Spitzer source 
matching, we request photometric quality flags for JK 



to be at least "B", i .e. we select only detections with 
SNR > 7 and photometric uncertainties below 0.1551 
mag. Relaxing this criterion by allowing selection of the 
flags "C" and "D" , we end up with the stellar disk frac- 
tion of 41/90 = (46 ±9)%. Among the BD population in 
our FOV, we find that the fraction of sources exhibiting 
MIR excess is 15/26 = (58 zt 19 ) %, in agreement with 



results fr om Jay awardh ana et al. (2003). The YSO sam- 
ple from Bontemps ct al.| ( |2U0l |""^includes both stars and 
BDs, down to ~ 0.05571^0. Combining our numbers for 
stars and BDs, we get the disk fraction of (45 ± 7)%, in 
agreement with |Bontemps et al.| ( |200T| . 

Similar results were obtained in other star forming re- 
gions. For a Ori, Chamaeleon-I, and IC 348, disk frac- 
tions are found to be ~ 40% and ~ 50 — 60% for low- 



mass s tars and brown dwarfs, respectively (Luhman et al. 
2008b| ). In NGC 1333, |Scholz "ciar('lJiT) find an upper 
limit on the BD disk fraction of 66%. a Ori, Cha-I, and 
IC 348 are somewhat older than p Oph and NGC 1333 
(2-3 Myr vs < 1 Myr), indicating that there is no signif- 
icant change in the disk fractions over the first 3 Myr of 
the cluster lifetime. 

6. CONCLUSIONS 

In this work we have built upon the large-scale optical 
and near-inf rared survey of the L1688 core in p Ophiuchi 
presented in Geers et al. (2011). From the optical-|-NIR 
photometry, 309 objects were selected as candidate sub- 
stellar cluster members. 58/309 of these objects were 
targeted for spectrosco py, and one was c onfirmed as a 
substellar mass object ( jGeers et al.| 2011). Here we re- 
port the second round of spectrosco]Dic follow-up, where 
we obtained spectra of additional 46 candidates from the 
same list. Among these, we find 4 candidate substellar 
objects with Toff < 3000K, 3 of which are for the first 
time reported here (SONYC-RhoOph-3, 6, 7). 
From MOIRCS, 2MASS and Spitzer photometry. 



Geers et al. (2011) established an additional sample of 
26 candidate substellar members with disk excess. 12 
of these appear as spectroscopically confirmed BDs with 
spectral type later than M5 in the literature. We ob- 
tained spectra of 17 Spitzer BD candidates, including 9 
previously confirmed BDs. We report 5 newly discovered 
very-low-mass probable members of p Oph, including 3 
BDs (SONYC-RhoOph-5, 8, 10). From the list of the 
previously confirmed objects, we re-confirm 7 out of 9 
with Teff < 3200K. 

In summary, we report the discovery of 8 new very-low- 
mass candidate members in p Oph with Toff between 2550 
and 3200 K. Six of these objects are classified as brown 
dwarfs, and comprise one fifth of the entire currently 
known substellar population in p Oph. While it might be 
premature to derive the substellar IMF in this cluster, as 
the high and variable extinction seriously hampers most 
attempts to detect objects at the lowest masses, we esti- 
mate that the upper limit on the ratio of the number of 
low-mass stars to the number of BDs is in line with other 
nearby star forming regions. By analysis of the MIR ex- 
cess among the spectroscopically confirmed members of 
p Oph, we find the disk fractions to be ~ 40% and ~ 60% 
for stars and BDs, respectively. This is again in line with 
the results from the other nearby star forming regions. 

The authors would like to thank the Subaru staff, es- 



SONYC - p Ophiuchi 
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Fig. 



8. — Spatial distribution of spectroscopically confirmed members in p Oph. Contours and lines are same as in FigurelTl with addition 
of the lowest extinction contour. Stellar members and BD are shown as black and red symbols, respectively. Crosses mark objects without 
disks, and diamonds those with disks (i.e. showing MIR excess). 
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